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0. Introduction 



Multiple scale techniques (MST) originated in Poincare works have been developed by 
many authors, mainly in solving (partial) differential equations related to physical prob- 
lems in celestial mechanics or in fluid dynamics. All these methods have a common 
mathematical purpose: to avoid resonances or secularities appearing in the usual or con- 
ventional perturbative theory. From a more physical point of view, one can see the MST 
as adaptable methods that feel the underlying physical phenomena in order to fit them. 
In other words, the usual perturbative theory tends to impose its choices while MST are 
flexible and compose with the real medium. 

In this work, we apply one of the various MST to the quantum anharmonic oscillator. 
Such studies have been initiated by Bender and Bettencourt (B&B) in two recent papers 
[|1], 0. They have found that the non resonance condition leads to a "mass renormal- 
isation" of the oscillator and - as a by-product - to the energy level differences of the 
quantum oscillator. This pioneering work was limited to the first non trivial order in 
MST perturbation of the coupling constant of the anharmonicity. The aim of the present 
paper is to extend this early study in several directions. First, we introduce an alterna- 
tive framework, which turns out to be more convenient than the B&B one for performing 
higher order calculations. Secondly, it turns out that we are able to obtain the energy 
levels themselves at these perturbative orders. In the third point, we show that the diag- 
onalization of the Hamiltonian is rather easy once the free Hamiltonian has been recast in 
an appropriate form. Finally, this approach leads to a natural and elegant method to find 
perturbatively the eigenvalues of the full Hamiltonian, far away from the original MST 
concept. 

The paper is organized as follows. In the first section, although the classical anhar- 
monic oscillator is studied in details in many textbooks ||, f§], we sketch some relevant 
points in order to further clarify the differences and the analogies between the classical 
and the quantum cases. In the second section we explain our framework and we work out 
the two first orders in MST perturbation, that includes the full solution of the Heisen- 
berg equations and the energy levels. The third section is devoted to general arguments 
showing that the method is compelled to work at any order, due to its connection with 
a certain unitary transformation which diagonalizes the Hamiltonian of the anharmonic 
oscillator. We postpone to an appendix some explicit results: the solutions of the Heisen- 
berg equations of motion and the energy levels of the full Hamiltonian, up to the order 6 
included. . 

1. The Classical case 

The classical anharmonic oscillator (CAO) is probably one of the most popular examples 
where the conventional and MST perturbative theories lead to obvious differences. Often, 
one speaks of Duffing equation instead, although this equation is nothing but the equation 
of motion of the CAO. To be precise, the Duffing equation is a second order non linear 
equation in the time variable, the solution of which being the position of the CAO. Starting 
from the CAO Lagrangian ( in units where the mass parameter is 1 ) 

L(q,q) = q 2 /2-u 2 q 2 /2-gq i 

one readily gets from the Euler-Lagrange equation: 



2 



(e) : q + u 2 q + Agq 3 = 0. 



The usual formal perturbation expansion reads 

oo 

q(t) = jT i g a q n {t), 

n=0 

( with some initial conditions, say q(0) = Q and q(0) = ), and the first equations one 
obtains from (e) are : 

(e ) : q + ^ 2 <?o = 0, 

(ei) : q\ + u; 2 gi = -iq^. 

Then the frequency of the solution of the homogeneous part of (ei) coincides with the 
frequency of go(^) = Qcosut, which generates a resonance in the solution qi(t) of the full 
equation (ei): 

Q 3 

qi(t) = — -(cos3u)t — cosut — 12ujtsinujt) . 

Hence qi(t) is unbounded and the truncated expansion q (t) + gq\{t) cannot be an 
acceptable approximation of q{t) for times t larger than u/Q 2 g, however small g may be. 
The flaw is even worse at the higher orders. It is obvious on this simple example that the 
perturbative solution develops spurious behaviour which is absent in the exact solution. 
Indeed, it is well known that the exact solution is bounded and periodic. 

The main idea of MST for dealing with this problem is the introduction of new vari- 
ables, independent and appropriate, and we refer to textbooks for an extensive review of 
the various possibilities. Here we concentrate on the anharmonic oscillator. Some meth- 
ods take into account ab initio that the circular functions play a major role. For instance, 
in the Poincare method, one looks for sine and cosine solutions whose argument is still ut 
but where u is now an arbitrary function of the coupling constant, actually lu = Y,g n u> n . 
Then one has to find the u^s , order by order, to discard the resonance. We do not 
insist on the application of these methods to the CAO because we believe they are not 
suitable for the quantum case. Another class of MST seems to be of a larger use, since 
there is no "prerequisite" in these methods. The MST we will use, also called Derivative 
Expansion Method, belongs to this class: it promotes the time variable to be a function 
of the coupling constant, namely t n = g n t. Actually, the method is not so rough and 
one first extends the function depending on t to an "extended" function depending on 
all the variables t n (n = 0, 1,2, ...) assumed to be independent ||. So, one introduces a 
new position function Q(T,g) depending on the collection T = {to, £1,^2, •••} of indepen- 
dent variables t n . This function is considered as an extension of the true position in the 
Lagrange formalism, which is recovered by restricting Q to the section t n = g n t of the 
T-space: q(t, g) = Q(T, g)\t n = g ^t- 

Then, forgetting temporarily any reference to the coupling constant in these t n vari- 
ables, one expands in power of g the new position function Q : 

00 

Q(T,g) = J2g n Qn(T). 

n=0 
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One obtains from (e) the following set of equations, limited here at the three first 
orders : 

D 2 Q (T)+u 2 Q (T)=0, 



DlQ x {T) + tfQ^T) = -2D D 1 Q (T) - 4Q^(T) , 



D 2 Q 2 (T) + uo 2 Q 2 {T) = -{D\ + 2D D 2 )Q (T) - 2D D 1 Q 1 (T) - l2Ql{T)Q x {T), 
using dj dt = Y, n g n D n , D n = dj dt n . 

The basic principle of the method now consists in adjusting the t\ dependence of Q{T) 
so as to eliminate the secularity in the second equation, next the t 2 dependence of Q(T) 
so as to eliminate the secularity in the third equation, and so on. We shall not work out 
the derivation here (it can be found for example in ref 3 or 4) and we merely give the 
solution up to the second order in g in its final form, for further classical versus quantum 
discussions : 

q(t, g) = |[exp(-i(m + &)) + ^(1 - HH) eX p(-3z(fit + b)) + ^ exp(-5z(m + &))] + C.C., 



n = co(l + ^-^),X= y -^ (1) 



where : 

3A 15A\ , ga 2 

), A = — 

2 16 ; u 

and a and b are two real integration constants fixed by the initial conditions (here un- 
specified). 

Since the (perturbative) energy is conserved, it can be computed most easily by choos- 
ing t = -b/tt or t = (| - b) /Q in q(t, g) : 

aV,, 9A 25A 2 . 
^ = ^-( 1 + T + -^ r ) + 0(A 3 ). (2) 

We conclude this first section by a few comments. As far as we know, all the multiple 
scale techniques dealing with the secularities of the classical anharmonic oscillator are 
successful. However this is not a general feature, and some methods are not suitable for 
certain problems. Moreover it is absolutely not our purpose to discuss on a rigorous basis 
the mathematical aspects of the secular or non secular perturbative expansions. 



2. The Quantum case : Derivation 

The quantum anharmonic oscillator (QAO) has been studied in the paper of B&B through 
the Heisenberg equation of motions for the relevant operators and we will follow this 
method. The main difference between the work of B&B and ours is that we will use the 
creation and annihilation operators to manage the problem of removing the secularities. 
At first sight the gain in doing this choice is not obvious and perhaps not essential. 
Moreover one can detect in the B&B paper an indication pointing to this direction. Let 
us look at the couple of equations (21) in their work [l|], which can be written as : 
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D X Y = -CX - XC and D X X = CY + YC 



where X and Y are ti dependent, self-adjoint operators while C is a constant, self-adjoint 
operator. The authors proceed with some arguments "suggesting" the form of the solu- 
tion, with the help of Weyl ordered products and Euler polynomials to deal with these 
equations. Of course, it seems difficult, or at least hazardous, to generalize at high orders 
a "suggestive" method, which could be seen as a reminiscence of the Poincare method, 
but we have more convincing arguments to leave this path. First, using Z = X + iY, the 
previous couple of equations reduces to the single equation : 

D 1 Z(t 1 ) = -i(Z(t 1 )C + CZ(t 1 )) 

whose solution is Z(ti) = exp(— iCti)Z(0) exp(— iCti), as it is easy to check. The op- 
erator Z(ti) is closely related to the creation/annihilation operators. Once derived the 
expression of the creation/annihilation operators, it is not necessary, in order to proceed 
further, to write down the position operator. Indeed, almost all the informations, the 
"mass renormalization" effect and the difference of energy levels, are already contained 
in the argument of the exponentials. Secondly, the Heisenberg equations in terms of cre- 
ation/annihilation operators are first order differential equations in place of the second 
order one for the position operator, which simplifies noticeabily the whole procedure. To 
be honest, one has the disadvantage to carry both creator and annihilator, but this is 
not a serious complication. Lastly, there appears also a large variation between the B&B 
works and ours in the status of the initial conditions: we do not use these conditions as 
in the classical case, which is the way taken by B&B. This point will become obvious 
throughout our study. 

We start with the QAO Hamiltonian H written in terms of the momentum p and 
position q operators in convenient units (H — cu — 1) : H = p 2 /2 + q 2 /2 + gq 4 , where g 
is assumed to be a "small" (positive) coupling constant. Whithin the Heisenberg picture, 
the dynamics is governed by the equations : 

q = i[H, q],p = i[H,p], 

supplemented by the canonical commutation relation [q,p] = i, valid at all times. The 
Heisenberg equations give : q = p and p = —q — 4gg 3 . Writing as usual q = (a + aJ)/y/2 
and p = —i(a — a^)/v^2, the Hamiltonian becomes: 

H(a,a\g) = 1/2 + a*a + g(a + a f ) 4 /4 (3) 

together with : 

[a(t,g), a(t,g)] = 1 , Vt, (4) 

where, to avoid possible confusion later on, we have kept track of the variables t and g. 
The Heisenberg equation for the annihilator : 

a(t, g) = i[H(a(t, g), a f (t, g),g), a(t, g)], 

reads, in our case : 
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a(t, g) = -i(a(t, g) + g(a(t, g) + a\t, g)) 3 ). (5) 
Since the Hamiltonian is conserved, its formal solution is: 

a(t, g) = exp(i#(a(0), a f (0), g)t)a(0) exp(-i#(a(0), a f (0), #)t), 
with a(0) = a(0, g). 

We now turn on the formal series of the multitime perturbative expansion, similar to 
that used in the classical case. First one introduces an operator valued function A(T, g) 
depending on the collection T of independent variables tj. This function is considered as 
an extension of the true annihilation operator in the Heisenberg picture, which is recovered 
through the restriction: 

a(t,g) = A(T,g)\ tj=gH . (6) 
Then the time derivative becomes : 

a(t,g) = Y,g n D n A(T,g)\ tj=gH . 



n>0 

Secondly, A(T, g) is expanded as : 

A(T,g) = J29 n A n (T). (7) 

n>0 

As for the initial conditions to be associated with the equation of motion (5), one 
notices that 

a(0,g) = Y,9 n A n (0). (8) 

n>0 

This forces us to choose between two possible starting viewpoints : 
either a) : a(0, g) is taken as independent of g , which implies 

A n (0) = 0,Vn> 1, (9) 

or b) : the previous condition is not imposed, in which case the initial values of a(t, g) 
must be considered as a function of g. 

It turns out that both approaches lead to consistent multitime expansions. In fact, 
the choice a) was (implicitely) adopted by B&B. However, these authors did not extend 
their analysis beyond the first order. In this paper, we rather follow the procedure b), 
which we found much more convenient, and in a sense, more natural. 

The equation of motion for a(t, g) gives us the following infinite system for the A n (T)'s 



n-l 



D A n + iA n = - D n- m Am - i QmQrQs (n = 0,l,2..) (10) 

771=0 m,r,s,>,0 

m+T — \-s — n— 1 

where Q n =A n + A' n , or explicitely : 

D o A o + iA o = 0, (lO.a) 
D A ± + iA ± = -DtAo - iQl (lO.b) 
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D A 2 + %A 2 = -(D 2 A + D 1 A 1 ) - t{QlQ l + Q0Q1Q0 + QiQ 2 ), (lO.c) 
etc... 

A simple check shows us that any formal solution of (10) generates via (6) and 
(7) a formal solution a(t,g) of (5). In particular, this implies that, for such a solu- 
tion, [A(T, g),A^(T,g)]\ t=g j t is independent oft. Of course, this does no mean yet that 
[A(T, g), A^(T, g)] is independent of T, allowing us to impose : 

{A(T,g),A^T,g)] = l,VT, (11) 

in order to insure the canonical commutation relation (4) . However, one can look for 
those solutions of (10) which are subjected to the stronger condition (11), if such so- 
lutions do exist indeed, i.e. if no inconsistencies or obstructions arise in their iterative 
construction. Together with (7), this entails : 

[A (T),Al(T)] = l 

VT (12) 

E n m=0 [Am(T),A^ m (T)] = 0, n>l 

We are now ready to construct step by step the resonance-free solution of the problem. 
To zeroth order, the equation (10. a) and the first equation (12) yield : 

A (T) = A i(T 1 )exp(-^ ) (13) 

with 

[A 01 (T 1 ),AUt 1 )] = 1,\JT 1 , (14) 

and the notation : Tk = tfc+i, ...}, (k = 1, 2, ...). 

Then, one can proceed to the first order step by inserting eq (13) into eq (10. b) : 

D A 1 + iA 1 = —{D 1 A 01 + iiA^A^ + A^iAn + 4i^oi)) exp(-it ) - i(A 3 01 exp(-3it ) 



+ 4i exp(+3it ) + (AllA 01 + 4iAn4i + A 01 All)exp(+it )). (15) 

Before integrating this equation, one has to get rid of the first resonant term on the 
right hand side, which would produce a contribution growing linearly with t (= t). This 
leads to the condition : 

D ± A 01 = -iiA^Al, + A 01 Al 1 A 01 + 4i^oi) (16) 

which will fix the t\ dependence of A Q1 . 

To do that, let us first introduce the self-adjoint operator N(T) = Aq(T)A (T). 
Thanks to (13) and its creator version, N(T) is only T\ dependent: N(T) = A^T^A^Ti). 
Moreover as a consequence of (14), A i(Ti)N(Ti) = (N(Ti) + l)A i(Ti). Lastly, from (16) 
, one observes that D 1 N(Ti) = 0. Thus is also independent of t\ and (16) can be now 
written in the tractable form : 

D\A m = -3iA)i(7i)JV(T 2 ), 

which produces : 

A 01 (T 1 ) = A 02 (T 2 ) exp(-3iiV(r 2 )i 1 ). 
7 



This allows us to write down the first order annihilation operator : 

A (T) = A 02 (T 2 )exp(-i(t + 3iV(T 2 )t 1 )). (17) 
At the same time, (14) becomes : 

[A) 2 (T 2 ),4 2 (T 2 )] = 1,VT 2 . (18) 

One can now come back to the form of (15) exempted of secularity to obtain its general 
solution : 

Ai(T) = Al^T,) exp(-3it )/2 - 4 3 1 (T 1 )e^(+3^ )/4 - 3iV(T 2 )4 1 (T 1 ) exp(+*t )/2 

+ C 1 (T 1 )exp(-it ). (19) 

where the operator Ci(Ti) is an integration "constant". The latter must be so adjusted, 
if possible, as to insure that the second equation (12) , 

[A (T),A\(T)} + [Ai(T) , Aq(T)] = 0, (20) 

be fulfilled at all times T. Here, it turns out that (20) is satisfied by taking simply 
Ci(Ti) = 0. One ends up with : 

Ai{T) = A 3 (T)/2 - Af{T)/A - 3iV(T 2 )4(T)/2 (21) 

and the first order step is complete. 

Before going further, some comments are in order. First, writing the position operator 
Qo + 9Qi, one notes that the coefficients of exp(±it ) in go get corrections coming from 
qi. It means, in the position formalism, the scheme used by B&B, that one would have 
to take into account the solutions of the homogeneous second order differential equation. 
Secondly, it appears in (21) that any power of exp(+it ) (resp. exp(— it )) is multiplied 
by the same power of Aq^Ti) (resp. A 01 (Ti)). Such a correspondance, which is specific to 
our way of managing the initial conditions, will be a guide throughout our study. Lastly 
the solution of the homogeneous equation in the classical case is different. This variation 
with the quantum case is due to the different status of the initial conditions. 

Clearly, one can go iteratively through the higher order steps by similar (although 
rapidly tedious) calculations as long as the integration "constants" analogous to C\(Ti) 
can be properly adjusted. As in the first order step, we gather in eq (10. c) the terms 
containing exp(—it ), since exp(— it ) is again (and always) solution of the homogeneous 
equation. Because D\Ai{T) does not provide such a term, we just have to take into 
account the non derivative part of the right hand side of eq (10. c). Through an intensive 
use of the relation A 01 (T 1 )N(T 2 ) = (N(T 2 ) + l)A 01 (Ti), this expression can be reduced 
to: — 3Ao 2 (T 2 )(17iV 2 (T 2 ) + 7)exp(— it — 3iN(T 2 )ti)/A, and the non resonance condition 
coming from the second order reads : 

D 2 A 02 (T 2 ) = 3tA 02 (T 2 )(17N 2 (T 2 ) + 7)/4. 

This equation shows that N(T 2 ) is in fact independent of t 2 , too, (i.e. N(T 2 ) = Al 3 (T 3 )A 03 (T 3 ) 
) and we find through integration: 
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A 02 (T 2 ) = A 03 (T 3 ) exp(+3*(17iV 2 (T3) + 7)t 2 /4), (22) 
whereas (18) becomes: 

[A m (T 3 ),4 3 (T 3 )] = l,VT !l . (23) 

Collecting equations (7), (17) and (22), we see that the non resonance conditions, up 
to the second order, imply that the first order term of the expansion of the annihilation 
operator is, in the variable t: 

a (t, g) = a (0, g)(exp(-tt(l + 3gN - 3g 2 (17N 2 + 7)/4)) + 0(g 3 )), (24) 

which exhibits a large difference with the classical case : 17 is a prime number, difficult 
to link with the other prime number 5 coming in the CAO frequency (1) . We will discuss 
later on this CAO/QAO (apparent) discrepancy. Nevertheless, the result, equation (24), 
is in perfect agreement with the perturbative expression of the energy levels of the QAO, 
as calculated by standard methods : 

E n (g) = 1/2 + n + 3g(l + 2n + 2n 2 )/4 - g 2 (l + 2ra)(21 + Yin + 17n 2 )/8 + 0(g 3 ). (25) 

Indeed, a straightforward argument based on the formal expression of a(t, g) in the 
Heisenberg picture shows us that the frequency appearing in (24) for N = n should 
coincide with E n (g) — E n _i(g). This is readily checked. 

Turning back on the second order equation (10. c) cleared from its resonant terms, we 
obtain its general solution : 

A 2 (T) = -15A 3 (N - l)/4 + 3^/16 + 3(23N 2 + 7)^/8 

+ 21(N- 1)4 3 /16 - 4V8 + C 2 (T 1 ) exp(-rto), (26) 

where A and N stand for A (T) and N(T 3 ). In contrast with Ci(Ti) in (19), the operator 
C 2 (Ti) cannot be taken as vanishing, because the second condition (12), 

[A (T),A\(T)} + [Ai{T), A\(T)\ + [A 2 (T) , Al(T)} = 0, 

would not be fulfilled. Imposing this and using eqs (12), (17), (18) and (26), one finds 
instead an appropriate expression for the solution of the homogeneous version of (10. c), 
namely : 

C 2 (T0 expHto) = -9A)(T)(1 - 3iV 2 )/32. (27) 

Let us notice that the (operator) coefficients of exp(±i£ ) which appear in the zeroth 
order solution get corrections from the first and second orders, and the coefficients of 
exp(±3it ) which appear at the first order get also corrections coming from the second 
order. Such a behaviour still holds at the third order, as we have checked. 

So far, the perturbative expression of the energy levels of the QAO (which was not 
our main goal) did not show up in full within our MST procedure. Yet, it can be found 
(without appealing to other perturbative methods) by inserting a(t,g) as given by equa- 
tions (6), (21), (26) and (27) in the Hamiltonian (3). Obviously, we are waiting for an 
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expansion in powers of g polynomially dependent on A and Aq, up to the second order 
in g : 

H = H + gH 1 + g 2 H 2 + O(g 3 ), 
The result is that the if/s are function of N — AqAq, not of A and Aq separately : 
H — 1/2 + iV + 3g(l + 2N + 2N 2 ) /4 - g 2 (l + 2JV)(21 + YIN + 17iV 2 ) /8 + 0(g 3 ). (28) 

This feature, which technically appears as an accident due to many cancellations, is in 
fact easy to understand. One observes, at each step, the to, h--- dependences of A (T) 
arise from the exponentials only. Since the t-dependence must eventually disappear from 
the conserved quantity H, a proper balance between A (T) and A (T) is expected in each 
of the monomials Hj, namely as many creators as annihilators. Then, whatever are the 
number and the order of the A 's and A 's in those polynomials, the commutatiom relation 
(12) allows us to cast the Hj's in the form of polynomials in N = A\Aq. Furthermore, 
anticipating a result to be proved in the next section, N(T) is not only independent of 
to, ti and ti but in fact of T altogether: N = Aq(0)A o (0). On account of the Heisenberg 
algebra (12) (taken at T = 0) this implies that the spectrum of N is the set of non negative 
integers, and the expression (28) for H is in complete agreement with (25) indeed. 

Finally, we have to explain the apparent discrepancy noticed earlier between the clas- 
sical and quantum results. Let us write the classical energy (2) for uj — 1: E c — 
a 2 /2 + 3ga 4 /8 + 0(g 2 ). Now for large n the quantum energy (25) reduces to : E q = 
n + 3gn 2 /2 + 0(g 2 ). Then the natural correspondance is a 2 / 2 — > n plus a quantum cor- 
rection so adjusted as to insure E c = E q + 0(g 2 ). One finds a 2 /2 — > n — 3gn 2 which, 
inserted in the classical frequency (1) gives Q = 1 + 3gn — 51g 2 n 2 /4, in agreement with 
the large n quantum frequency from (24), derived within the MST scheme. 

3. The Quantum case : General discussion 

In the previous section, in particular on eqs (21), (26) and (27), one observes that the 
construction is made with two elementary bricks A (T) and Aq(T), where A (T) is the 
first term in the MST expansion of the annihilation operator. We have also pointed out 
the simple connection between the operator N and the Hamiltonian. More precisely, the 
first perturbative results exhibit the following features : 

1) N = aI(T)A (T) is independent of t , t u t 2 ..., i.e. of T. 

2) The "nonhomogeneous" parts of A n (T) depend on T through the basic operators 
A (T) and Aq(T). The same is true for the "homogeneous" parts C n {Ti) exp(— it ) which, 
after determination of C n (Ti), can be recast in the form of functions of A Q (T) and A (T) 
only. 

3) The operators H and N commute. 

If these features persist at all orders, then (putting aside any consideration of conver- 
gence) one should obtain in the limit : 

oo 

A(T,g) = Aq(T) + J2g n F n (MT)A(T)) (29) 

n=l 



10 



where : 

oo 

A (T) = A expH(t + E fn(N)t n )) (30) 

n=l 

together with [Aq, Aq] = 1, and where the F n 's are some polynomial functions of A and 
A\ while the /„'s are some polynomial functions of N. We will show below that the res- 
onance -free solutions of the perturbative multitime equations of motion do exist indeed, 
and have the general form (29)- (30). This means, in particular, that no obstructions 
are encountered in determining the integration "constants" C n (Ti) and giving them the 
appropriate form. Eqs (29) and (30) then yield : 



a(t, g) = a (t, g) + Yl 9 n F n {a {t, g), c&(t, g)) (31) 

n=l 

where : 

oo 

a (t, g) = A exp(-i(l + ]T g n f n (N))t) (32) 

n=l 

and N = AqA is a constant operator. 

Actually, these facts result from the full equivalence between the iterative process 
described in the previous section and the perturbative determination of an unitary trans- 
formation which brings the Hamiltonian to a diagonal form. 

In order to prove this equivalence, let us consider the spectral decomposition of the 
Hamiltonian : 

oo 

H{a Q ,a\,g) = 1/2 + ala + g(a + 4) 4 /4 = E n(g)\n, g >< n,g\, 

n=0 

where {| n, g >} is the orthonormal basis made of the "perturbed" eigenvalues of H 
(for future convenience, a is written here as do). We also introduce the "unperturbed", 
orthonormal Fock basis {| n >} induced by the operators a and %, together with the 
unitary transformation which maps the former onto the latter : 

\n>= U(g)\n,g> (n = 0,1,2...). 

The unitary operator U(g) is determined up to a N dependent, arbitrary, right phase 
factor, where iV = ajao. 

Then if we define H d as H d (a ,a , g) = U\g)H(a,o,a\,g)U(g), we have : 

H d (a , o{, g) = H(a(g), a\g),g), (33) 
where we denote by a(g) the annihilation operator transformed by U(g) : 

a(g) = U\g)a,U{g). (34) 

We also have : 

oo 

H d (a ,al,g) = J2 E n(g)\n >< n\. 

n=0 
11 



In other words, the unitary transformations (34) of the dynamical variables is that 
one which diagonalizes the Hamiltonian in the {\n >} basis. The perturbative form of 
eqs (33) and (34) are 



a(g) = «o + 9 n a n (35) 

71=1 

and 

OO oo 

H d (a , al g) = Y 9 k H k = - + a\a G + £ g k H k , (36) 

n=0 1 n=l 

where the explicit expressions of the H k s in terms of the obtained by substituing 

the perturbative form (35) in H (a(g) , a) (g) , g) , and expanding. For the QAO Hamilto- 
nian we are interested in, thoseiffc's are: 

k 

H k = Y a ln a k-m + Y q m q r q s qi/4:,(k = l,2,....), 

71=0 m,r,s ,£>Q 

m + r + s + £ — k — 1 

where q m = a m + a^. 

The operators a n (n = 1,2,3...) in (35) are determined recursively as polynomial 
functions of ao and by requiring that : 

i) U(g) be unitary indeed, or equivalently (due to Von Neumann theorem ||) that the 
commutator of the "new" variables a(g) and a\g) be canonical: [a(g),afl(g)] = l,V<7,i.e.: 



Y Wm, a n -m] = $n,o, (ra = 0, 1, 2...). (37) 

m=0 

ii) Hd be diagonal indeed in the {| n >} basis or, equivalently, that [Ha,N] = 0,\/g, 
i.e. : 

[H k ,N} = 0,(k = 1,2,3..) (38) 

Eqs (38) implies that all the H k s are functions of the operator N = ala only. In 
particular these operators H k commute between themselves. 

Evidently, the equations (37) and (38) must admit solution for {a n }, due to the mere 
existence of the unitary operator U(g) and its formal perturbative expansion. However 
there is no uniqueness property because of the phase freedom in the mapping U(g). In 
our case of QAO with standard quartic interaction, the "minimal" solution {a n } is such 
that each a n is a polynomial of degree 2n + 1 in a and a\ with rational coefficients, and 
monomials of odd degrees only : 

«n= J2 (x nAl a l Q N k + y nAl N k al l ). 

1 < l cdd< 2 "+ 1 

Let us now define the T dependent operators a n (T) by: 

oo oo 

a n {T) = exp(i ^ H k t k )a n exp(-z ^ H k t k ), (n = 0, 1, 2, ...) (39) 

fc=0 fc=0 
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where t k = g k t. We claim that these operators obey exactly the canonical commutation 
relations (12) and the differential equations (10) which serve previously to determine the 

Because of (37) , this is immediate for the relations (12). As for the equations (10) , 
one first derives from (39) : 

oo oo 

D n _ m a m {T) = i exp(i ^ H k t k )[H n _ m ({a r }), a m ] exp(-i H k t k) 

k=0 k=0 

which, summing up, yields 

n n 

(T) = i ^ [H n _ m ({a r (T)}),a m (T)\, (40) 

m=0 m=0 

where the commutativity of the H k s has been used twice. On the other hand, using (7) 
together with (33) and (36) to express H(A(T, g), A^(T, g),g) in terms of the function H k , 
one readily finds that the equations (10) read as well: 

n n 

]T D n _ m A m (T) = i ]T [tf n _ m ({A r (T)}), A m (T)}, (41) 

m=0 m=0 

identical to (40). This is actually true not only for the QAO but for a general interaction. 

Therefore, {A n (T)} can be identified as one of the solutions {a n (T)}, which establishes 
the equivalence of the two schemes, and hence the consistence of the multitime method 
we used, together with the general validity of the assertions 1) to 3) put forward at the 
begining of this section. 

It is possible now to comment on the "mass renormalisation" introduced by B&B. Since 
H d is a pure function of N, H d = H(N), eq(39) for a (T)\ tj=gjt = A (T)\ tj=gH = a (t,g) 
reads: 

a Q (t,g) = exp(+i#(A0t)a o exp(-i#(iV)t), 
or, by using clqN — (N + l)ao : 

ao(t, g) = a exp(-i(H(N) - H(N - l))t). 

Together with (28), this gives: 

ao(t, g) = a Q exp(-it(l + 3gN - 3g 2 (7 + 17iV 2 )/4) + 0(g 3 j). 

The "renormalisation" phenomenon can be pinned down to the fact that H(N) — 
H(N — 1) is the trivial identity operator at the zeroth order, and becomes a true operator 
for higher orders. 

As mentioned at the begining, and apparent on eqs (29) and (30), the solution A(T,g) 
constructed there corresponds to initial conditions depending on A and g . If one insists 
in having the perturbative solution with prescribed g— independent initial condition : 
A(0,g) = a, with [a, a^] = 1 , this is easily achieved by a few additional manipulations. 
Indeed, it is sufficient to invert order by order the relation : 

oo 

a = A + J2g n Fn(A ,4), 

71=1 
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(which is straightforward in spite of the non commutative algebra) to get 



A = a+ Y,9 n G n {a,a^), 

n=l 

and to reinsert this expression for A in eqs (29) and (30), as well as in N = AqA , 
truncated at the relevant order. Then, of course, the expression of A(T,g) in terms of a 
and has no longer the "simple" structure that it exhibits in terms of A and Aq. 

To conclude this section, we wish to stress again that the arguments presented there 
are quite general, not specific of the QAO. If one considers an Hamiltonian which is 
the sum of an harmonic oscillator one and a "potential" represented by a self-adjoint 
operator function of the position and the momentum, such an analysis can be repeated. 
Actually, the equivalence between MST and unitary transformation diagonalizing the 
Hamiltonian is likely to be a rather general feature. In particular, the previous discussion 
can be extended in a rather straigthforward way to systems with more than one degree 
of freedom. 

Furthermore, the equivalence between the multitime approach and the perturbative 
construction of the relevant unitary transformation must have a classical counterpart. In 
the classical framework, multitime expansions should appear as essentially equivalent to 
the construction of appropriate canonical transformations, following the Poincare - Von 
Zeipel method |7j], or some of its disguises. As a matter of fact, one can find indication of 
such a connection in the literature |4], ||. This aspect of the question, which we have not 
touched upon in the present paper, might deserve a further study. 

4. Conclusion 

In this paper, we have used the anharmonic oscillator in the Heisenberg picture as a 
model for investigating the practicability of the Derivative Expansion Method, of com- 
mon use in classical physics, within the quantum framework. This method turns out 
to be successful in providing us with the perturbative expansion of the time dependent 
dynamical variables together with the energy levels, which we have derived explicitely up 
to the second order. We also have proved that this MST is equivalent to the perturbative 
construction of an unitary transformation diagonalizing the full Hamiltonian, leading to 
a step-by-step algorithm for the calculation of the previous quantities at any order, and 
thereby strengthening the status of the Multiple Scale Techniques in quantum mechanics. 
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Appendix 



We give below, up to order 6: 

i) the coefficients a n of the expansion (35) of the annihilator a(g) in terms of N — ajao, 

ii) the coefficients Ek(n) of the expansion of the energy levels: 

1 oo 

E n (g) = - + n + Y,9 k Ek(n)- 

Z n=l 

Both have been computed by applying the algorithm described in section 3 (eqs (37) 
and (38)). 

i) To be simpler, a and stand for ao and aj. 

ai = (2a 3 - at 3 - ; 

a 2 = (-9a+120a 3 +6a 5 -120a 3 A^+27aA^ 2 +84at-42a t3 -4a t5 +42ATa t3 +276Ar 2 at)/32 ; 

a 3 = (6092a 3 + 756a 5 + 8a 7 - 8844a 3 iV + 4422a 3 iV 2 - 378a 5 iV + 1278aiV - 1062aiV 3 
-1708a t3 - 464at 5 - 6a t7 - 9282iVa t + 2406iVa t3 + 232iVa t5 - 1203iV 2 a t3 
-9042iW)/128 ; 

a 4 = (-200645a + 1546416a 3 + 380868a 5 + 9264a 7 + 40a 9 - 2975280a 3 iV + 2143296a 3 iV 2 
-714432a 3 iV 3 - 322896a 5 iV + 80724a 5 iV 2 - 3088a 7 iV - 500298aiV 2 + 162755aiV 4 
+506760at - 358200a t3 - 221832at 5 - 6696at 7 - 32at 9 + 673392iVa t3 + 186464iVa t5 
+2232iVat 7 +3040992A^ 2 at-472788iV 2 at 3 -46616iV 2 at 5 +157596iV 3 at 3 +1365240A^ 4 a t )/2048; 

a 5 = (116798776a 3 + 51228696a 5 + 2189520a 7 + 20832a 9 + 48a 11 - 266946576a 3 iV 
+255315936a 3 iV 2 -121842648a 3 iV 3 +30460662a 3 iV 4 -58614828a 5 iV+24750360a 5 iV 2 
-4125060a 5 A^ 3 -1300128a 7 iV+216688a 7 A^ 2 -5208a 9 A^+52602092aA^+41073824aA^ 3 
-6417388aiV 5 - 21539684at 3 - 28787584at 5 - 1542096at 7 - 16320a t9 
-40a tn - 12162525(Wa t + 50282976iVa t3 + 32551232iVa t5 + 912600iVa t7 
+4080iVat 9 - 47389884iV 2 at 3 - 13618080iV 2 a t5 - 152100iW 7 - 2199146767VV 
+22248396iVV 3 + 2269680iV 3 a t5 - 5562099iV 4 a t3 - 55675938iV 5 at)/8192 ; 

a 6 = (-2649077789a + 19854323040a 3 + 14799326898a 5 + 1000498176a 7 + 15874840a 9 

+78720a 11 + 112a 13 -15744a 11 iV-52410470592a 3 A^+59605775856a 3 iV 2 -37821182832a 3 A^ 3 
+13464443 160a 3 iV 4 - 2692888632a 3 A^ 5 - 20808622800a 5 N + 11852140500a 5 iV 2 
-3324992400a 5 N 3 +415624050a 5 N 4 -817924896a 7 iV+242212752a 7 iV 2 -26912528a 7 iV 3 
-7253184a 9 iV + 906648a 9 iV 2 - 16271788323aiV 2 - 609787599 laiV 4 + 521267535aiV 6 
+4255953324at - 2581523304a t3 - 8101045372at 5 - 691648560at 7 - 12242240a t9 
-64720at n - 96at 13 + 12944iVa tn + 7571823000iVa t3 + 11245609120iVat 5 
+562441728iVat 7 + 5584128iVat 9 + 35458238196iW - 9129805056iW 3 
-6326409960iV 2 a t5 - 165946104iV 2 at 7 - 698016iV 2 a t9 + 5783860872iV 3 a t3 
+ 1757503840iV 3 a t5 + 18438456iW 7 + 29695249 188iW - 2055444390iV 4 a t3 
-219687980iV 4 a t5 + 411088878iV 5 a t3 + 47684835487V 6 a t )/65536. 

ii) 

E 1 (n) = 3(1 + 2n + 2n 2 )/4 ; 

E 2 (n) = -(1 + 2n)(21 + 17n + 17n 2 )/8 ; 
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E 3 (n) = 3(111 + 347n + 472n 2 + 250n 3 + 125n 4 )/16 ; 

E 4 (n) = -(1 + 2n)(30885 + 49927n + 60616n 2 + 21378n 3 + 10689n 4 )/128 ; 

E 5 (n) = 3(305577 + 1189893™ + 2060462n 2 + 1857870n 3 + 1220765n 4 + 350196n 5 
+ 116732n 6 )/256 ; 

E e (n) = -(1 + 2n) (65518401 + 146338895n + 213172430n 2 + 139931868n 3 + 85627929n 4 
+18794394n 5 + 6264798ri 6 )/1024. 

Several number theoretic properties of the Ek(n) 7 s are worth pointing out. First, 
all the coefficients Ck p of n p in Ek{n) are rational and positive, and the signs of the 
Ek(n) J s alternate, as it should be. Perhaps new are the following observations: whereas 
the denominator in the expression of Ek(n) is a power of 2, the numerator is always a 
multiple of 3 (for integer n). This peculiarity was already noticed by Bender and Wu || 
for the ground state (n = 0). It thus turns out to hold for the excited levels too. Also, 
the sum of the numerators of the coefficients Ck p in each Ek(n) is a multiple of 5. Finally, 
if one expresses the Efc(n)'s in terms of the variable m = n + |, one observes that they 
are even polynomials with positive coefficients (multiplied by — m if k is even). More 
than that, all the zeroes of these polynomials are pure imaginary. This means that all the 
zeroes of Ek(n) lie on the line n = — | + iy ! 
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